Theory of non- Fermi liquid and pairing in electron-doped cuprates 
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We apply the spin-fermion model to study the normal state and pairing instability in electron- 
doped cuprates near the antiferromagnetic QCP. Peculiar frequency dependencies of the normal 
state properties are shown to emerge from the self-consistent equations on the fermionic and bosonic 
self-energies, and are in agreement with experimentally observed ones. We argue that the pairing 
instability is in the d^2_y2 channel, as in hole-doped cuprates, but theoretical Tc is much lower than 
in the hole-doped case. For the same hopping integrals and the interaction strength as in hole-doped 
materials, we obtain Tc ~ lOK at the end point of the antiferromagnetic phase. We argue that a 
strong reduction of Tc in electron-doped cuprates compared to hole-doped ones is due to critical 
role of the Fermi surface curvature for electron-doped materials. The d2,2_j,2 -pairing gap A(k, oj) is 
strongly non-monotonic along the Fermi surface. The position of the gap maxima, however, does 
not coincide with the hot spots, as the non-monotonic d^2_y2 gap persists even at doping when the 
hot spots merge on the Brillouin zone diagonals. 

PACS: 74.25.-q, 74.20.Mn 



I. INTRODUCTION 

The fascinating properties of high-temperature super- 
conductors continue to attract high interest of condensed- 
matter community over the last two decades. Most 
of the extensive experimental and theoretical studies of 
the cuprates have been performed on hole-doped mate- 
rials, such as La2_2;Sr2;Cu04, YBa2Cu306+a; (YBCO), 
Bi2Sr2CaCu208+:r(Bi-2212), Tl2Ba2Cu06-H:r, etc. In re- 
cent years, however, there has been growing interest 
in the properties of electron-doped cuprates such as 
Nd2-a;Ce^Cu04 and Prz-a^Ce^^CuO^. 

Hole-doped and electron-doped cuprates are in many 
respects similar. For both, optical conductivity mea- 
surements at small doping show a charge-transfer gap of 
about 2eVA4, and the electronic Fermi surface, measured 
by ARPES in electron-doped cuprates is reasonably de- 
scribed the same combination of hopping integrals as in 
hole-doped cuprate^. This impHes that both electron- 
doped and hole-doped materials are Hkely described by 
the same underlying Hubbard-type model^. The normal 
state behavior of optical conductivity near optimal dop- 
ing is also very similar in hole-doped and electron-doped 
materials ^*i. Most importantly, the superconductivity 
in both types of materials has d^2_y2 symmetry, as evi- 
denced by, e.g., ARPES measurements of the momentum 
dependence of the pairing gapS*^. 

On the other hand, the phase diagrams of hole-doped 
and electron-doped cuprates are somewhat different. 
Electron-doped cuprates have a much wider range of an- 
tiferromagnetism than hole-doped materialsiSiii. There 
is a strong evidence that the Mott physics is not very rel- 
evant in electron-doped cuprates near optimal doping as 
the optical data show^ that the 1.7eV charge-transfer gap 
completely melts away as the electron doping approaches 
its optimal value x ~ 0.15. There is also no analog 
in electron-doped cuprates of the pseudogap phase be- 
tween the antiferromagnetic and superconducting phases 



on the phase diagram. Some pseudogap behavior in the 
normal state has been detected in Ndi.85Ceo.i5Cu04, but 
the onset temperature for this behavior clearly tracks the 
Neel temperature ~, and so a (narrow) pseudogap phase 
is likely just a magnetic fluctuation regime of a quasi- 
2D antiferromagnetic. In the absence of the pseudogap, 
the phase diagram of electron-doped cuprates resembles 
a "typical" quantum-critical phase diagramJA - the super- 
conducting Tc{x) forms a dome above the antiferromag- 
netic quantum-critical point (QCP)^^. 

The superconducting properties of the two classes of 
materials also differ. First, Tc in electron-doped cuprates 
is around 10 — 20K which is almost an order of mag- 
nitude smaller than in hole-doped cuprates, and up to 
a factor 500 smaller than the pseudogap T* , which, as 
a large group of researchers believe, is the onset tem- 
perature of the pairing without coherence in hole-doped 
cuprates. Second, the pairing gap in electron-doped 
cuprates varies non-monotonically along the Fermi sur- 
face - it increases at deviations from Brillouin zone di- 
agonals, passes through a maximum, and then decreases. 
The non-monotonic behavior of the d— wave gap has been 
originally introduced to explain Raman data2i. Later, the 
non-monotonic gap has been directly observed in ARPES 
experiments^. In hole-doped cuprates, the rf— wave gap 
is monotonic, within error barA 

The large difference between Tc in electron and hole- 
doped cuprates (and even larger difference with T* for 
the hole-doped materials) despite apparently similar hop- 
ping integrals and the strength of the Hubbard interac- 
tion calls for an explanation. In this paper, we argue 
that the primary difference between the pairing instabil- 
ity temperature in electron-doped and hole-doped ma- 
terials is that in hole-doped cuprates the pairing pre- 
dominantly involves antinodal fermions, while for elec- 
tron doped cuprates the pairing involves fermions near 
zone diagonals. We show that for the pairing of near- 
diagonal fermions, the Fermi surface curvature has a very 



strong and negative effect on Tc, and reduces it by more 
than two orders of magnitude. This eventually leads to 
Tc ~ lOK in electron-doped cuprates. 

To solve the pairing problem, one needs to select a 
pairing mechanism. Following earlier works by one of 
us and the others ^ffi^Ti^fJ ^ we assume that the pairing is 
of electron rather than phonon origin, and that a dom- 
inant pairing interaction between fermions is mediated 
by collective, Landau-overdamped magnetic fluctuations. 
The idea of spin-fluctuation pairing was earlier applied 
to hole-doped cuprates, with the assumption that Mott 
physics does not play a substantial role in the normal 
(not pseudogap) state, and becomes relevant only well 
inside the pseudogap phase. For hole-doped cuprates, 
the Luttinger Fermi surface is such that the hot spots 
(Fermi surface points separated by the antiferromagnetic 
momentum) are located near the corners of the Brillouin 
zone. The pairing near a QCP predominantly involves 
fermions from around the hot spots, and the region of 
pairing instability forms a dome on top of the antiferro- 
magnetic QCBi&. The instability temperature increases 
with underdoping and saturates at T* ^ 0.0391&, where 
g is the effective interaction, which is comparable to the 
charge transfer gap. Using 1.5eV for g yields T* ~ 500K, 
which is consistent with the onset temperature for the 
pseudogap behavior. The relation between T* and the 
actual Tc is a separate issue which we will not discuss in 
this paper. 

We will apply the same itinerant spin-fermion model 
to study the spin-mediated pairing in electron-doped 
cuprates. Like we said, the electron-doped cuprates seem 
even more suitable for an itinerant description than the 
hole-doped materials as there is no non-magnetic pseu- 
dogap, and the charge-transfer gap melts away near opti- 
mal electron doping. We assume, following RPA studies 
of the static spin susceptibility, that antiferromagnetism 
sets in at around the doping when the hot spots merge on 
the Brillouin zone diagonalsiii^fliSi. At this doping the 
antiferromagnetic Brillouin zone touches the Fermi sur- 
face at the diagonal points ki? = (7r/2, 7r/2) (see Fig. ^, 
so that 2ki? coincides with the antiferromagnetic wave 
vector Q. At smaller dopings, the system is magneti- 
cally ordered, and the fermionic Fermi surface displays 
electron and hole pockets. 

We first consider normal state properties near a 2kp 
antiferromagnetic QCP. We reproduce and extend earlier 
result of Altshuler, loffe, and MillisiS2, that at QCP, the 
self-energy of a nodal fermion has a non- Fermi liquid form 
and scales as S(ti;) ex w". The exponent a is close to unity 
though itself varies with frequency. In this situation, the 
normal state behavior is close to that in a marginal Fermi 
liquid. 

We next consider the pairing of these non-Fermi liq- 
uid fermions and show that at QCP, the system is still 
unstable towards (1^2 _y2 gap opening, i.e., there is a 
dome of the pairing instability around critical doping 
where antiferromagnetism sets in. The pairing instabil- 
ity temperature at QCP (which for electron-doped sys- 



tems we will label as Tc because of the absence of the 
pseudogap) still scales with g but because of the strong 
de-pairing effect associated with the Fermi surface cur- 
vature, Tc ~ O.OOO65 for the actual curvature of the 
cuprates Fermi surface. This is 50 times smaller than 
0.03g for hole-doped cuprates. The uncertainty of Tc 
due to the approximate nature of the estimate of the 
curvature is actually very weak as Tc as a function of the 
curvature is almost flat in a wide range of the curvatures. 

We also show that the pairing gap is non-monotonic 
along the Fermi surface, and passes through a maximum 
at some deviation from the zone diagonal. The location 
of the gap maxima does not track hot spots as at 2kF 
QCP, the hot spots are right along the zone diagonals. 

That an antiferromagnetically mediated d^2_y2 pairing 
survives in a situation when the hot spots are along the 
zone diagonals is not intuitively obvious, for the strongest 
pairing interaction involves quasiparticles for which the 
dj.2_y2 superconducting gap vanishes. However, one can 
easily see that the separation between the near-nodal 
fermions with opposite signs of the (1^2 -y2 gap is on aver- 
age closer to Q than that between fermions with the same 
sign of the gap (see Fig. fflb). Because of this difference, 
there is still attraction in the da,2_j^2— wave channel for 
the spin-fluctuation mediated interaction. Furthermore, 
we will see that for the pairing at QCP, the kernel of 
the gap equation actually does not have any extra small- 
ness associated with the hot spots location along the zone 
diagonals, i.e., Tc is formally of order g. This happens 
because the gap equation relates A(k) near the two di- 
agonal hot spots ]i.F and hp -I- Q on the Fermi surface, 
where both gaps are linear in the deviation from the di- 
agonal. Then the gap equation effectively becomes an 
equation on the slope of A(k) and the latter one does 
not contain any extra smallness, as we will see. 

The dj.2_y2 pairing near 2kp instability has been ear- 
lier considered within BCS theorj*2i, where the authors 
found that d— wave Tc remains small but flnite at QCP. 
Our results agree with Keim^ in that the d— wave attrac- 
tion survives at QCP, but we argue that the pairing in- 
volves non- Fermi liquid fermions, and that Tc scales with 
the upper boundary of the quantum-critical regime (mod- 
ulo small prefactor) rather than with the upper bound- 
ary of the Fermi liquid behavior. The strong reduction 
of Tc in electron-doped cuprates compared to the onset 
temperature for the pairing in hole-doped cuprates was 
also obtained in FLEX studies of magnetically-mediated 
pairingii (although the difference was less drastic than 
in our case). We believe that the numerical results of 
RefJi describe the same physics as our analysis. Non- 
monotonic variation of the d— wave gap also agrees with 
recent numerical studies^i. 

A short summary of our results was presented earlier 
inSi. 

The paper is organized as follows. In Section ^J we 
briefly review the spin-fermion approach. In Section ITTTI 
we study the properties of electron-doped materials in 
the normal state. In subsection IIII Al we develop pertur- 



bation approach to fermionic and bosonic self-energies. 
Self-consistent solution of the Dyson's equations on the 
fermionic and bosonic self-energies is covered in subsec- 
tion ^^B| In subsection IIII (^1 we use the obtained self- 
energy to find the frequency dependences of the conduc- 
tivity and Raman absorption. In Section ITvl we address 
the pairing problem and obtain Tj, and the momentum 
dependence of the gap. In IIV Al we consider the limit 
when the parameter representing curvature of the Fermi 
surface in the gap equation formally tends to zero. In 
lIVBI we consider the case of finite curvature. We present 
numerical solution in IIV B l1 and an approximate analyt- 
ical solution in IIV B 2l and find that the strong variation 
of Tc with the curvatures can be well understood. In Sec- 
tion we sum up the main conclusions. Some technical 
details are relegated to Appendices. 



II. MODEL 

In this paper we apply the spin-fermion model to a 
particular form of electron dispersion in two dimensions. 
In general, low energy dynamics of strongly-interacting 
electrons may be studied by integrating out the high- 
energy part of electron-electron interaction. Close to a 
collective instability, the corresponding bosonic collective 
mode tends to become gapless, and has to be included 
in the low-energy model. The interaction with the near- 
gapless collective mode alters the fermionic dynamics and 
may lead to non-Fermi-liquid behavior. The same inter- 
action also gives rise to the pairing. 

A spin-fermion model has been previously applied to 
hole-doped cuprates, and we refer the reader to earlier 
literaturei^ for the details and the justification of the 
model. Here we start right with the Hamiltonian: 

k q 

+3Xl'=k+q,a'^"/3Ck/3Sq, (1) 

k,q 

Here Ck are the low-energy fermions and Sq are their col- 
lective bosonic fiuctuations in the spin channel. The spin- 
fermion interaction is described by the effective coupHng 
g\ self-consistency of the model requires g to be smaller, 
or, at most, of the order of the bandwidth g < W. 
Static spin susceptibility xo(q) is determined by the high- 
energy physics. We assume that it has Ornstein-Zernike 
form near the antiferromagnetic vector Q —{tt^tt): 



Xo(q) = Xa/iC 



Qf) 



(2) 



The correlation length ^ parameterizes proximity to the 
instability. At QCP which we only study below, ^~^ = 0. 
The interaction g and xo appear only in combination 
9 — 9^X0 which is the actual (measurable) spin-fermion 
interaction. 

While the static spin susceptibility is an input, the low- 
frequency dynamic part of the susceptibility (01 should be 




Figure 1: Left: Fermi surface at the antiferromagnetic QCP 
with 2kp = (7r,7r). Diamond-shaped dashed lines bound the 
magnetic Brillouin zone. The diagonal points of the Fermi 
surface (nodal points of the d^2_y2-wave gap) now become 
"hot". Right: Graphic explanation of the attraction in the 
d^2_y2 channel: parts of the Fermi surface on the same side 
of the zone diagonals are on average closer to the Q separation 
than the parts on the opposite sides, leading to attraction in 
the d^2_y2 channel (plus and minus are the signs of the d^2_y2 
gap). 



determined within the model, together with the fermionic 
self-energy. The reasoning is that the damping of the 
bosonic fiuctuations occurs only by a decay into particle- 
hole pairs with energies smaller than the bosonic fre- 
quency. Because ^ is peaked at the antiferromagnetic 
wave vector Q, the interaction between fermions and 
bosons will predominantly involve fermions near the hot 
spots (Fermi surface points separated by Q). 

The computational tool to solve the model of Eq. (QJ 
at large S, is Eliashberg-type theory in which one solves 
for fully renormalized fermionic and bosonic self-energies 
while simultaneously neglecting vertex corrections and 
the self-energy dependence on momentum transverse to 
the Fermi surface. This computational procedure is sim- 
ilar, but not equivalent to FLEX. Eliashberg theory at 
strong coupling has been discussed extensively in early 
literatureiS and we will just use it in our work without 
further discussion. 



III. NORMAL- STATE PROPERTIES 
A. Perturbation approach 

The system behavior at 2kF QCP within Eliashberg 
theory has been studied by Altshuler et ali^. They used 
RG procedure to sum up the most divergent diagrams 
and obtained fermionic and bosonic self-energies at the 
lowest frequencies in the normal state. We obtained very 
similar results using a somewhat different computational 
procedure (see below). We also obtained fermionic self- 
energy at intermediate frequencies which are mostly rel- 
evant for the pairing problem. 

A set of self-consistent Dyson's equations for the 
fermionic and bosonic self-energies in the normal state 



IS 



X-\q,u;) = xoiq)-' +Ii{q,u;), (3) 

G-\k,uj) = Go(k,cj)-i+iS(k,cj). (4) 

The Green's function of noninteracting fermions is 

GQ^(k,ui) =iuj- eu, (5) 

and ek is the electron dispersion. Bosonic self-energy 
(polarization operator) contains contributions from non- 
umklapp and umklapp scattering between hot spots sep- 
arated by Q. Hereafter we denote by q the deviation of q 
from the antiferromagnetic wave vector Q unless other- 
wise noted. In these units, the total polarization operator 
is 2(n(q, ft) + n(-q, 17)), where 

n(q, n) = 2g^ J ^^G(k, c.)G(k + q,u; + n). (6) 



Then 



x(q,^) 



Xo 



q2 + 2xo(n(q,r!)+n(-q,r!))- 



(7) 



The peculiarity of a spin-fermion theory near a 2kF in- 
stability is that the Fermi velocities at the hot spots are 
antiparallel, and to keep n(q, il) finite, one has to ex- 
pand the spectrum of the fermions in the vicinity of these 
points up to the second order in the momentum compo- 
nent along the Fermi surface: 



fk+Q 



-vpk^ + P'^kl 



(8) 
(9) 



Here the momentum components kx, ky are normal and 
tangential to the Fermi line respectively, and k is mea- 
sured from the hot spot (see Fig. QJ; /? parameterizes the 
curvature of the Fermi line, k = 2/3^/wf- The radius of 
curvature used in2^ is k^, — 1/k. The special case of a 
nested Fermi surface has been analyzed in -S A. Virosztek 
and J. Ruvalds, Phys. Rev. B 42, 4064 (1990). 

For noninteracting fermions with the Green's function 
©, Eq. © yields 



no(q,f7) = 



2'kvfI3 



VL^ 



El 



E. 



q' 



where 



Eq = -VFQx + f3^ql/2. 



(10) 



(11) 



Already this zeroth-order expression Ijl0|l differs in two 
important ways from previously studied case^^ when the 
two hot spots separated by Q had nearly orthogonal ve- 
locities. First, at q = Q Il{il) scales as V" instead of the 
conventional \fl\, and, moreover, diverges when /? — *■ , 
so curvature of the Fermi line cannot be neglected. Later 
we will see that /? 7^ is even more important for the 
pairing. It is interesting to trace how ljl()|l converts into 



a conventional Landau damping when the hot spots are 
moving away from the zone diagonals. This is done in 

Appendix El 

Second, a conventional Landau damping \n\ term rep- 
resents an anomaly in analytical properties of bosonic 
self-energy. External frequency fi appears as a factor 
because it regularizes the integral over internal frequen- 
cies by separating the poles into opposite semi-planes. 
At the same time, the polarization operator (I10|l gives 
a fractional power of external frequency not as a result 
of anomaly, but simply because ft sets an effective cut- 
off for the ultraviolet-divergent integral over internal uj. 
This distinction becomes relevant at a finite T: while 
Landau damping is a purely quantum phenomenon, and 
so does not depend on temperature, the self-energy lfT?Hl 
possesses u>/T -scaling. Evaluating the polarization oper- 
ator for free fermions but at a finite T, we obtain 



no(q,fi) 



9' 


f°° v/V4w2 + f^2 + 2tj 


2ttvfP J 


-00 4Tcosh „r^' 



= VTf{n/T,Eq/T), 



duj{U) 
(13) 



where 



fia,b) 



9 



vV4p + P + 2a; 



da:. 



2ttvf(3 7-_oo 4cosh^((a; - b)/2) 
At q = 0, 17 = 0, the polarization operator is simply: 



(14) 



n, 



— 



2ttvfI3 



|C(i)|(^/2-l) = const VT, 



(15) 



where C{^) ~ —1.46 is the Riemann Zeta functional. 
When T ^ 0, the function l/4Tcosh^(a:/2r) in the de- 
nominator of ltT2|l becomes a (5(a;)-function, and ltT2|l con- 
verts to ifTHji . 

The electron self-energy is given by 

^^^x(k-k',c.-c.')G(k',c.'), 

(16) 
The coefficient 3 = (tct in II16|I comes from the three com- 
ponents of the bosonic spin fiuctuations. In Eliashberg 
theory, the self-energy depends on frequency and on the 
momentum component along the Fermi surface. We pa- 
rameterize the position of the quasiparticle on the Fermi 
surface by the momentum component ky (see Fig. QJ, i.e., 
label S](k, w) as S(fcy,a;) along the Fermi surface. This 
does not imply that kx is zero as k^ and ky are related: 
from Ek = 0, kx = —(5'^ky/vF along the Fermi surface. 
Since fermions are fast compared to bosons, the integral 
over momenta in IJ16II can be factorized: the one perpen- 
dicular to the Fermi surface involves only fast electrons, 
and the one along the Fermi surface (ek' = 0) involves 
slow bosons. The corrections from keeping the momen- 
tum transverse to the Fermi surface in the bosonic prop- 
agator is of the same smallness as the vertex correction 
(in this, Eliashberg theory differs from FLEX where the 
momentum integration is not factorized) 



Substituting 



Gik,Lu)-^ = iiu + E{k,Lu)) - ek 



(17) 



into itTfijI and integrating over momenta transverse to the 
Fermi surface, we obtain for the self-energy 

I](fc„c^)=y^ — J dk'yxik,k'y,n), (18) 



The reduced bosonic propagator 



x{ky,k'y,n) 



■iirvp 



X(k'-k,f])|,^^,^,^0 



(19) 



is taken between the two points k and k' on opposite 
sheets of the Fermi surface. 

The reduced propagator depends on the curvature in 
two ways. First, the polarization operator 11 depends 
on curvature via the overall factor and the En term in 



Eq. JTTJ. Second, there is a direct /3 dependence in the 
static part of the susceptibility, as {kx — k'^Y term in 
X reduces to /?^(fc^ + k'y'^)'^/vp once we use the Fermi 
surface relations k^ = —f}'^ky/vF and k'^ — P'^k'y^ /vp- 

For the self-energy exactly at the hot spot (7r/2,7r/2), 
the first dependence is much more relevant than the sec- 
ond one. Indeed, setting the external momentum fc = 
and taking k' to be at the Fermi surface we find that Ec^ 
and i?_q terms with q = k' reduce to E]^, = -fi'^k'y^ /2, 
and -E-k' = 3/3^A:^^/2, i.e., i?±k' scales with the curva- 
ture. The polarization operators n(q, J7) and n(— q, Vt) 
given by Ill0|l depend on the ratio ri/i?±k', and there- 
fore the dependence of i?±k' on the curvature is relevant. 
On the other hand, the P'^k'y'^/vp term in Xo^ i^ clearly 
subleading to k'J^ at small ky. 

Keeping the dependence on the curvature only in the 
polarization operator we obtain that away from the QCP, 
when ^~^ 7^ 0, S(fcy = 0,w) = 2](w) has a usual 
Fermi liquid expansion in powers oi uj. At the QCP 
the Fermi liquid behavior breaks down, and S acquires 
non-Fermi-liquid frequency dependence. To find it, we 
assume and then verify that typical ky in the integral 
for the self-energy are much larger than typical fl, and 
expand the denominator in Q in frequency. Using 
Eq = -P'^k'y^/2 and E.^ = ZP^k'y^ /2 valid when k' is 
at the Fermi surface, we find that the expansion holds in 
k'y" + V3<Zo|A:;| + |r!|gg/4wo|fc;i, where 



90 = g/T^VF, 

UJQ = {(j(3/2TTVFf'- 



(20) 

(21) 



The extra static term -v/Sgol^yl coming from 11 just re- 
fiects the fact that once k' is nonzero, the actual distance 
between the diagonal Fermi surface point and a Fermi 
surface point on the opposite sheet of the Fermi surface 
is not exactly (tTjTt), so the susceptibility acquires an ad- 
ditional "mass" term. Note also that the leading term 
in the frequency expansion has a conventional |ri|/|/Cy| 
term, typical for small momentum scattering. 



Substituting this expansion into x(9 = k'^) and in- 
tegrating over k'y in ifTSJI , we obtain that the momentum 
integral is logarithmic and is cut from below by frequency 
Q. (this justifies the assumption that typical momenta are 
larger than Vl). Integrating finally over frequency, we ob- 
tain 



E(^) 



-— wlog — , 

47r luq 



(22) 



This marginal Fermi Hquid, cjlogo;, behavior of the self- 
energy at small frequencies was first detected in2^. A 
simple analysis shows that it extends to frequencies of 
order wq. At larger frequencies lo > loq, typical k'y in the 

integral for the self-energy are of order go(^/^o)^^^- At 
these momenta, E±q ^ (fiwo)^/^ < Q, so the polarization 
operator can well be approximated by its zero momentum 
form n(0, fl) (X vO. Substituting this form into Ijl8|l . we 
find 



S(c.) 



1/4| |3/4 

ujJ \uj\ ' sgnuj. 



(23) 



At finite ky, the form of the self-energy T,{ky,uj) is 
rather involved. Below we will only need self-energy at 
Lu > ujq since in electron-doped cuprates Wq ^ lOmeV is 
small (see below). Parametrically, the leading depen- 
dence comes from the polarization operator, via vf<1x 
term in Eq, which for the particles at the Fermi surface 
reduces to 0^ {k"^ + k'y^) ~ (i'^ky as {ky-k'yY is small and 
irrelevant. Substituting this into n((7, 0) and evaluating 
the self-energy, we obtain 



Y^{ky,{jj) 



,l/4|,„|3/4 



sgnix) 



1 



p-'-kl 



(24) 



where a = Oil). This formula shows that the self- 
energy begins decreasing at deviations from the nodal 
direction ky ^ ^/ 0^- At the same time, the momentum 
dependence in 1I2H1 is very weak from practical point of 
view. The /cj,"* term in the static part of x~^ leads to a 
stronger 1/(1 + [ky/koYY''^ dependence, but this hap- 
pens at larger ky as fco oc uj^'^. We assumed that the 
self-energy scaling behavior is roughly in between the two 
analytical dependences - the self-energy Y,{ky,uo) rather 
weakly depends on ky up to some ky which approximately 
scales as w^'^, and smoothly decreases at larger ky. 

Finally, for the calculations of Tc in Section IIVL we 
will need the self-energy at a finite T. We evaluated it 
numerically by replacing integration over frequency by 
summation over discrete Matsubara frequencies. 

i;(fcy,w,„) = sgnWmT ^ dkyX{ky,k'y,Un), (25) 

|o„|<w„ 

where x is the same as in H19|l . Numerically we found 
that the explicit temperature dependence of ^{ky,ujm) is 
again weak at w > wq, so the self-energy at finite temper- 
ature has to a good accuracy the same functional form 
as at r = 0, but in discrete Matsubara frequencies. 



B. Self-consistent self-energies 

The marginal Fermi liquid form of the self-energy at 
Lu < ujo and the ui^''^ dependence at oj > uiq were ob- 
tained using the bare fermionic propagators. We need 
to verify whether the results survive when calculations 
of the fermionic and bosonic self-energies are performed 
self-consistently, using the full propagators. 

Landau damping linear-in-fJ term in the polarization 
operator is an anomaly, and it does not depend on 
whether the polarization bubble is evaluated with bare 
or full fermionic propagators, as long as the self-energy 
depends only on frequency. However, as we already 
discussed, the fl^^^ term in 11 is not an anomaly, and 
the form of n(q = 0, 51) does generally depend on the 
fermionic self-energy, which in turn depends on the func- 
tional form of n. This implies that lowest-order results 
are not sufHcient, and one has to carry out full self- 
consistent calculations. This was first noticed in—. 

Fortunately, these calculations are not necessary at 
uj > LOo as at these frequencies the self-energy E ex 
^(wo/w)^'^ is smaller than uj, and the renormalized 
fermionic Green's function is close to the bare one. In 
this situation, the expression n(q, Q) ex fi^/^ that was 
obtained using free fermions is a good approximation for 
interacting fermions as well, and no corrections are nec- 
essary. Next, using the fact that the self-energy S(a;) is 
expressed via the fermionic density of states, and the lat- 
ter does not depend on the fermionic self-energy (again 
as long as S depends only on w), we find that Eq. Ij23|l 
survives. Vertex corrections and other corrections be- 
yond Eliashberg theory can be treated in the same way 
as for other quantum-critical problemsiSi^SiS^, and are 
irrelevant. 

On the other hand, the marginal Fermi liquid form 
of the self-energy at frequencies smaller than ujq, Eq. 
1I22I1 . does not survive in the higher-order corrections. 
In cuprates ujq is small, and so this is not very rele- 
vant to our study, but it is interesting to address it 
from the principle point of view and we briefiy discuss 
it. Altshuler et al" found that next order term yields 
S2(w) ex wlog^w and conjectured that higher-order log- 
arithmic corrections form geometric series that exponen- 
tiates to a power law 



S(^) 



cjq|w|^ "sgnw. 



(26) 



They found a = \/3/47r « 0.14. We obtained a nearly 
identical result using a different computational proce- 
dure with does not require that the series is geomet- 
ric. Namely, we obtained a self-consistent solution for 
the self-energy iflfijl by evaluating the polarization opera- 
tor n with the full Green;s function Ijl7|l . We assumed 
that S(a;) ex |a;|^~", where a is unknown, re-evaluated 
the polarization bubble with the full fermionic G [this 
gives n(q,0) ex £;^^~"^/^", ll{0,n) ex f](2-a)/2]^ g^^sti- 
tuted the result into the integral for S(u;), and solved the 
self-consistent equation on a. The calculations, which are 



detailed in Appendix El yield a « 0.15 which is surpris- 
ingly close (although not identical) to a = 0.14 obtained 
by Altshuler et &&. 



C. Conductivity and Raman Response 

We now use the result for the self-energy to compute 
optical conductivity and Raman response in the normal 
state. Complex conductivity (j{iLj) can be computed from 
the Kubo formula: 



CTapii^n) 



tU) 



1 



tU) 



0), 



^^^^(^c.„)-^l^^(^0) 



(27) 



where Ilrj is the current-current correlator with zero 
transmitted momentum. Dependence on real to is 
found by the transformation iw„ -^ co + iO from 
Matsubara frequencies cj„. Assuming constant veloc- 
ities vf, the current-current correlator Ir^UiUn) = 
e^VpSaf3Tl{iLOn)/2g'^ is proportional to the Matsubara 
bubble n(iw„) with zero transmitted momentum. Ne- 
glecting corrections to the current vertices (which is jus- 
tified when S predominantly depends on to), we obtain 

anu!) ex / dky / — r — -. 

^m J Jo UJjn + 2:[ky,UJ,n-e)+^{ky,e) 

(28) 
At zero temperature and finite frequencies which we con- 
sider, the short-circuiting of the conductivity considered 
in^ is not important, and the conductivity can be ap- 
proximated by expanding the denominator in II28|I in 
the self-energy, integrating the self-energy over ky, and 
putting the result back into the denominator. The largest 
contribution to conductivity then comes from the nodal 
regions where the self-energy has a non- Fermi liquid form. 
To integrate the self-energy explicitly over ky requires 
substantial computational efforts as there are two distinct 
sources for ky dependence in E (see previous subsection). 
We carry out approximate calculations - we assume that 
S(fcy, e) is independent on ky and equal to ^{ky = 0, e) = 
S(e) for \ky\ < ky^-^^^ and then falls off rapidly, such that: 



Cf{iuj) ex 



'-i/jinax 



(Wm) 



W„ 



de 



Y.{uJr, 



WY 



(29) 

Following the discussion in preceeding subsection, we 
assume that for uJm > 1^0 , the threshold momentum 



^y,niax 



(Wm) 



1/4 



We also assume by continuity that 



for LU < Wo, fc!y,max(Wm) ~ ^f, 



Since ky 



depends 



on frequency, there is an ambiguity whether we should 
choose the external uj or internal e for the cutoff. One 
also should be careful when transforming from Matsub- 
ara frequencies to the real ones. We comment on this 
matter in more detail in Appendix El The conclusion 
is basically that this ambiguity is irrelevant for the fre- 
quency dependence of conductivity. 

At very small frequencies, w < wq, the self-energy dom- 
inates fermionic dynamics S((jj) > lo, and both 3cr and 




(o/eo„ 



Figure 2: Normal state conductivity as a function of frequency 
shows a scaling behavior for luq < u; < 40aJo. Left: |o"(tj)|, the 



behavior is indistinguishable from the cr{uj) 



depen- 



dence. Right: arctan3a(cij)/Kcr((i;), an almost constant value 
means that both Qa{uj) and Kcr(a;) scale as lj""'®*. 



^a have the same power-law behavior: <j{uj) cx lu' 



-1+2q 



-0.7 



. At larger frequencies S(a;) < lu, and one should 
not generally expect dia and Sct to scale with each other. 
Surprisingly, the scaling behavior extends, with almost 
the same exponent - we found that <t{uj) oc w~°-^'' over 
a very wide range up to ~ 40wo (see Fig. [21 . Such 
power-law behavior is not indicative of quantum-critical 
scaling, but rather a consequence of the flattening of 
the fermionic self-energy at high frequencies^, i.e. that 
5](a;)/ci; oc (wq/w)^''^ is a slow decaying function. The 
io~'' behavior of conductivity with 7 ~ 0.68 has been 
observed in Pri.85Ceo.i5Cu04 below 400meV ^. Both 
the exponent and the experimental frequency range are 
quite consistent with our results. A very similar behav- 
ior of the infrared conductivity at intermediate energies, 
also caused by the flattening of S, has been discussed for 
hole-doped cuprates^. 

Raman absorption measures the imaginary part of 
the fully renormalized particle-hole susceptibility at van- 
ishingly small incoming momentum, weighted with Ra- 
man form factors that depend on the scattering geom- 
etry. In the Big geometry the Raman form factors are 
jBigiM) OC (cosfca; — cosfcj,). They are the largest for 
fermionic momenta near (0, tt) and symmetry related 
points. In the i?2g geometry, the Raman form factors 
are 7^23 (k) oc sin kx sin ky , and the Raman signal comes 
from around {tt/2, 7r/2) points which are close to the diag- 
onal hot spots on the Fermi surface. For electron-doped 
cuprates these are the regions where we found non-Fermi 
liquid behavior of the fermionic self-energy. 

Unlike the current vertex for conductivity, the Raman 
vertex is renormalized by the interaction even when S((jj) 
depends only on frequency. The renormalization of the 
Big Raman vertex is not that relevant because the in- 
teraction peaked at or near Q — {it, tt) connects regions 
of the Brillouin zone where the bare vertex changes sign. 
The ladder series pertaining to vertex renormaHzation 
is then alternate in sign and roughly renormalizes jBig 
to 7_Bi /(I + A-jBin), where A is positive. On the other 



hand, the B2g Raman vertex jB2g (k) does not change 
sign under k — > k-l-Q. To a good approximation, we can 
approximate this vertex by a constant 753^ . The renor- 
malization of the Raman vertex then coincides with that 
of the density vertex, and like a density vertex, is then 
related to the self-energy by the Ward identity: the full 
7^2'' (^) ~ 7^2g (l+9cjS(ci;)). The vertex then leads to the 
additional frequency dependence 7B2g(^) ^ d^'S{uj) such 
that 7B2g(^) fx cj^" for lu < cvo and lu^^/'^ for uj > ujq. 

Convoluting two Raman vertices with the polarization 
bubble we obtain, approximately 



RB2g{^) ^ fcy,max(^) (1 + d^^{uj)) 



LU 



OJ - 



S(w) 



(30) 



Substituting the forms of the self-energy and ky^^s^^{Lu) r^ 
iu" for LU < uiQ, ky^i-aaxi^j) ^ Lo^ ^ ^ for tjj > LUo, we find that 
at T = 0, i?s2g(<^) is flat: it is a constant at a; < luq, and 
slowly crosses over to the (w/wq)"^ behavior at w > luq. 
The almost flat form of the Raman intensity at fre- 
quencies w > Wo is consistent with the experimental 
data2^. The experimental behavior at small frequencies 
is not flat, but according to2^ it is dominated by temper- 
ature effects which we do not consider here. 



IV. PAIRING PROBLEM 

The Hnearized gap equation at QCP is obtained using 
Eliashberg technique for collective-mode mediated pair- 
ing. Assuming that the pairing occurs in the singlet chan- 



nel, we write the pairing vertex as ^a/sip) 
and obtain (see Fig. \^ 



'a/3 



^(p) 



$(A;) = -3g2 ^ $(fc')G(fc')G(-fc')x(fc - k') (31) 



where k stands for a 3-component vector (k, cj), 
Ep = TY,^^(Pk/{2TTf, the Green's functions are full: 
G^^(fc) = Gq^ + S(fc), and the prefactor —3 comes from 
the convolutions of the Pauli matrices a^ia'Craa'CgB' — 



a'f3' 



— 3cr^^. The pairing vertex <l>(fc) is related to pairing gap 
A(fc) as $(k,w) = A{k, uj)Lu/{Lu + E{k,uj)). The solution 
of this linearized gap equation yields both Tc and the mo- 
mentum and frequency dependence of $(k, cj). At least 
close to Tc, this momentum dependence must be close to 
that of the true pairing gap. 

As in the normal state analysis, we use the fact 
that bosons are slower than fermions and approximate 
x{k — k') by Fermi surface to Fermi surface interaction, 
i.e, put both k and k' on the Fermi surface. Like be- 
fore, we also assume that the fermionic self-energy does 
not depend on the momentum component transverse to 
the Fermi surface. We will, however, keep the momen- 
tum dependence of S along the Fermi surface. We will 
see below that T^ very strongly and non-analytically de- 
pends on the Fermi surface curvature k = 20^ jvp. This 



strong non-analytic dependence comes from the depen- 
dence on curvature of the Fermi surface to Fermi surface 
interaction x(fc — k') (see below). Accordingly, we will 
keep the dependence on the curvature in x(^ ~ k') and 
in the fermionic self-energy (because S is an integral of 
x{k — k')), but neglect it in the Jacobian for the transfor- 
mation from the integration over (Pk to the integration 
over dck, i.e., set dkxdky = {l/vF)dekdky (the notations 
are the same as in Fig. 1). Under these assumptions, the 
momentum integration transverse to the Fermi surface 
involves the product of two fermionic propagators and 
reduces to 



j G{p)G{~p)de^ = 



\uj + Y.{ky,u})\ 



(32) 



Substituting this into Ij31|l and using the fact that the 
interaction involves momentum transfers near Q, we ob- 
tain 



$0(fca,Wn)--T 



where the kernel 



dk'yKn 



i{ky,k')<^Q{k' ujm), 



J^nml'^i 1^ I — 



X(fc,fc^CJn -UJm) 



(33) 



(34) 



The reduced x is defined in ((Tflll . ky and k'y are small, and 
the notations <I>o and $q imply that the pairing vertices 
on the right and left hand sides of Ij33|l belong to the 
opposite hot spots. 

Because the kernel lIMjl is invariant under the sign in- 
versions of frequency and momenta, the solutions are 
symmetric or antisymmetric in both. A solution with 
dx2_y2 symmetry is antisymmetric in ky. It is also sym- 
metric in frequency as at Tc the largest eigenvalue of 
the kernel reaches unity, and from a known theorem 
the eigenfunction of the largest eigenvalue is symmet- 
ric in a*2i. As a further simplification, we consider a 
"antiferromagnetic" dx^-j,^ order parameter made out of 
n = 1, n = 3, etc. partial components of the Big 
representation of the D^h group for tetragonal symme- 
try (i.e., coskx — cos ky, cosSk^ — cos 3ky, ...). For 
the odd n solutions, the pairing vertex (and the gap) 
change sign when the momentum is shifted by Q, i.e., 
^o{ky,ujn) — —^Q{ky,LL)n). Note that a "conventional" 
coskx — cos ky solution falls into this category. Using 
this condition and dropping the subscript we get from 



^{ky,LUn)=T}_^ / dk'yKnrn(ky,k'y)<^{k'y,UJ.m). (35) 

The remaining input is the polarization operator which 
is a part of x- We verified a posteriori that the dominant 
contribution to the pairing comes from uj > luq, where 
n((7, SI) scales as Q^'^ and weakly depends on q. In prin- 
ciple, one should use the full, finite T form of of the polar- 
ization bubble. However, the explicit temperature depen- 
dence of n only complicates the calculations but does not 
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Figure 4: Pairing instability temperature Tc in units of 
the coupling constant 'g vs. the curvature parameter r = 
g[3 /-KVp. Solid line gives the result of the numeric solu- 
tion, which requires increasingly long computational time 
as r — > 0. Dashed line is the analytic small-r form Tc = 
0.95gr/{l + 5r^^'^f. 



lead to any new physics. We assume without further dis- 
cussion that the T dependence of the polarization bubble 
does not alter Tc in any substantial way and will us e the 
same functional form Il{Qm) = (5^/27rwi?/3)-\/|f^m| as at 
T = 0, but use discrete Matsubara frequencies at a finite 
T. For accuracy, in numerical calculations we kept the 
full T dependence of the fermionic self-energy. We found, 
however, that the full form of 'Et{ky,ujm,T) differs very 
little from the zero-temperature form S(fcj,, tUm) taken at 
discrete Matsubara frequencies 

We stick to TV = 2 in this section. It is convenient to 
measure frequency and temperature in units of wq I|21|I as 
the fermionic self-energy at fc^ = is ll23ll . and measure 
momenta ky in units of go dOt- Switching to these new 
units, we reproduce Eqs. I|33ll and IjliJII with 



xik,k',n) = 



1 



2N {k ~ k'f + r^{k^ + A:'2)2 + ^\ 



The dimensionless quantity 



r = gP'^/nvj 



( 9 



\ 2'KVppp 



npF 



(36) 



(37) 



is proportional to both the curvature and the interaction. 
Both fc2 and fc"* terms in x(fc, k' , J7) come from the original 
(k — k')2 term in x(k — k') when we use that for the 



two fermions on the opposite sheets of the Fermi surface 
kx = P'^ky/vp and k'^ = — (3'^ k'y^ / v p ■ The normaUzation 
scale 0^0 also scales with r: 



"°=4^'' 



(38) 



and hence the curvature is present both in the over- 
all factor, and in the pairing susceptibility (and in the 
fermionic self-energy at a non-zero fc^, by virtue of the 
r— dependence of the susceptibility). 

The presence of the curvature both in the overall nor- 
malization factor for T and in the pairing susceptibil- 
ity is what qualitatively distinguishes the d— wave pair- 
ing of near-nodal fermions from the pairing of antinodal 
fermions, which was previously studied in the context of 
pairing in hole-doped cuprates. In the latter case, the 
pairing instability temperature is finite already at zero 
curvature, and the effect of the curvature on T^ is likely 
small (although this issue has not been analyzed in de- 
tail yet). In the present case, Tc just vanishes without 
curvature, so the curvature, parameterized by r, plays 
the major role. The same r also specifies momentum 
dependence of $(fcj,,ix'„), and, hence, of the pairing gap. 

Below we will solve the gap equation separately for the 
cases r — !■ 0, when only the r— dependence of the overall 
factor matters, and for finite r, when the rfc^ term in the 
pairing susceptibility has to be taken into consideration. 
We will see that the r— dependence in the denominator 
of l(3fijl becomes relevant already at very small r, and for 
r ^ 0.1, relevant to the cuprates (see below), the actual 
Tc is much smaller than the one obtained by keeping the 
curvature only in the overall factor. 



A. Solution at r ^ (vanishing curvature) 

When we put r = in the susceptibility li;ffi|l . the gap 
equation IJ35II simplifies considerably. First, the fermionic 
self-energy loses its dependence on ky and becomes 



Yiyky^ujm) — S(lj„j) 



sgnLj„i-7rT 



E 

n„i<|[j 



\^r. 



^1/4 



.1 



(39) 

We used the full summation in the numerical calculations 
below, but this is actually not even necessary - we found 
that for T = Tc, the functional form of the self-energy is 
close to the zero-temperature expression sgnwmlf^mP^'* 
for all Matsubara frequencies. 

Second, without the r term in x, momentum de- 
pendence of the kernel lIMIl is of a convolution type: 
K{ky,k'y) = K{ky - k'y). This implies that Eq. ^ is 
local in real space coordinate y conjugate to the momen- 
tum ky, as can be easily seen by taking reverse Fourier 
transform in ky. The kernel K{ky — k'y) is not double 
integrable: 



\K {k, k')\ dkdk' = oo 



(40) 



and so the integral equation in ky is non-Fredholmian, 
which means it does not have a countable spectrum with 
integrable eigenfunctions. However, the kernel is still uni- 
formly finite 



\K{k,k')\dk 



\K {k, k')\ dk' <oo, (41) 



and thus may have finite solutions. Two such solutions 
can be easily guessed: a symmetric one: $(fcy,a;„) — 
<i>(cj„) and an antisymmetric one: ^{ky,ijjn) = ky<^(ujn)- 
As functions of real space coordinate y these two solu- 
tions correspond to (5(j/)$(a;„) and —{dS{y)/dy)^{uJn), 
respectively, where S(y) is a delta function. The d— wave 
solution we are interested is the one antisymmetric in ky. 
In both cases, integrating explicitly over momentum and 
substituting the result for the susceptibility, we obtain 
one-dimensional integral equation for $(w„): 



$(c^„) = -ttT 



E 



$(Wm) 



k„-^m|l/4(|l 



(4 



Approximating the self-energy by the cu^'^ form simplifies 
^ to 



$(c.„) = -7rT^. 



<^(LUr, 



S(w„ 



11/4 ■ 



(43) 



This equation falls into a generic class of local Eliashberg- 
type gap equations with the effective local pairing inter- 
action xi{^) ^ '^~^ for which the fermionic self-energy 
scales as uj^~^ (i2i, see also Appendix EJ- Eq. (|l3ll cor- 
responds to 7 = 1/4. 

Note that: (i) the gap equation is fully universal and 
parameter-free (we recall that T is in units of wq), and 
(ii) that self-energy term is not small and cannot be ne- 
glected. At small LOm = 7rT(2m -I- 1) < 1, the kernel of 
the gap equation scales as l/w. This power is a combina- 
tion of l/w^/^ from the effective local pairing interaction 
and o;'^/'* from the fermionic self-energy. At LOm > 1, the 
kernel decays as l/w^'*, i.e., faster than 1/w, and the 
frequency sum converges even without '^(ujm)- Ref»2i ar- 
gued that, this gap equation has a non-zero solution at 
T = 0(1). 

We solved this equation numerically (the details are 
presented in Appendix and found, in actual units 



Tc{r ^ 0) w 6tJo 



R 9 



(44) 



The scaHng with the upper cutoff of the quantum-critical 
behavior (wq in our case) is the same as in hole-doped ma- 
terials (where the pairing involves predominantly antin- 
odal fermions), however, we emphasize that in our case 
the scale wq is by itself proportional to the Fermi surface 
curvature, while the corresponding scale for antinodal 
pairing remains finite in the absence of the curvature. 
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Note in passing that the prefactor is much larger in our 
case than for antinodal pairing, where T « 0.2 in units 
of the upper cutoff of the quantum-critical behavior i&. 



B. Solution at finite r 



1. Numerical solution 



Although Tc given by 14411 is finite, the solution of the 
gap equation without the r'^ky term in the susceptibil- 
ity is unphysical since the pairing gap continuously in- 
creases with kyi $(fcj,,CL'„) ex ky. This behavior refiects 
translational invariance of the pairing potential (or, in 
mathematical language, the non-Fredholmian property 
of the integral equation). When the r^k'^ term is kept in 
xiky, ky,uj), Eq. Ij^-{5|l becomes Fredholmian, its solutions 
integrable, and the gap vanishes at large ky. A numer- 
ical solution of the gap equation can then be obtained 
by standard technique. We solved Eq. ^l^^ numerically 
in the quadrant uj,k > 0, searching for a solution $„(A:) 
that is symmetric in w„ and antisymmetric in k. The 
semiinfinite integral over k was mapped onto the interval 
[0, 1) using the transformation x — (1 — fc)/fc, and then in- 
tegration was approximated by Gauss-Legendre quadra- 
ture: /p f{x)dx = '^if{xi)wi, where Xi and w^are the 
abscissas and weights of the quadrature. This procedure 
is called the Nystrom method^^. The integral equation 
then reduces to an (infinite in the Matsubara frequency 
index n) set of algebraic equations: 

This set was truncated and solved by standard LAPACK 
routines. This procedure gives good approximation to 
the largest eigenvalues of the kernel, and for the critical 
temperature we need just the largest one. 

We found that the effect of the curvature on Tc is very 
strong: above r > 0.001, the actual Tc rapidly becomes 
much smaller than lEjI . We plot Tc{r) in Fig. \^'m units 
of g/47r. We see that over a wide range 0.01 < r < 0.1, 
Tc ~ 0.00055 forms a plateau and weakly depends on r. 
Using the same g ~ 1.6eV and the t-t' dispersion as in 
hole-doped materiala^^ but with positive chemical poten- 
tial to reproduce the Fermi surface in Fig. 1, we obtained 
r ~ 0.08, which is within the region where Tc is almost 
constant, and Tc ~ lOK. For hole-doped cuprates, for 
the same parameters, the onset of the pairing was earlier 
estimated at Tc ~ O.Olg ~ 200-KiS,, although this num- 
ber may also be reduced by the curvature of the Fermi 
surface. 

The theoretical value of Tc ~ lOK at a magnetic 
QCP in electron-doped cuprates is in reasonable agree- 
ment with the experimenfeiiii^iS^. Maximum Tc in 
Nd2-2:Cea;Cu04 and Pr2-a;Cea:Cu04 is about 20- 25K, 
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Figure 5: The pairing vertex ^{ky,iJn) for r = 0.05 vs 
uj„/(2-kTc) and ky/qo {qo is defined in the text). Observe 
that ^{ky,uj„) is non-monotonic in ky. 



but, according tOi^^, Tc initially increases once the sys- 
tem becomes magnetically ordered. In any event, even 
50% agreement is quite reasonable given the number of 
approximations in our theoretical analysis. Note that if 
we used Eq. (|iill for Tc instead of the correct result, 
we would have obtained a very large Tc ~ 600K for the 
same parameters. This shows that the Fermi surface cur- 
vature truly plays a major role when the pairing involves 
near-nodal fermions. 

The reduction of Tc in electron-doped cuprates, com- 
pared to hole-doped cuprates with the same interaction 
strength was also obtained in FLEX calculationsii, al- 
though the reported difference was less drastic than in 
our analysis. 

In Fig we present the result for the momentum de- 
pendence of the pairing vertex at various frequencies. We 
see that the pairing gap is a non-monotonic function of 
kyi it initially increases with ky, but then passes through 
a maximum and decreases at larger ky . This result agrees 
with earlier BCS calculations^^, but we emphasize that 
our gap equation is of non-BCS form. We also emphasize 
that the position of the maximum in ^{ky,uj) is discon- 
nected from the location of the hot spots which in our 
QCP analysis are exactly along the Brillouin zone diag- 



onals, i.e., at k,. 



0. 



2. Toy model 

The large discrepancy between Tc{r) and Ij44|l can also 
be understood analytically, by expanding Tc in r beyond 
the 0{r) term. This expansion is rather non-trivial as 
one has to expand around a solution which diverges at 
large ky, and the divergence must be cut by the cor- 
rections which then obviously are not small at large ky. 
Since this effect is unrelated to the summation over Mat- 
subara frequencies, which at Tc is confined to first few 



m, we can simplify the model by dropping the frequency 
summation, and instead of II35|I . analyze the approximate 
form of the gap equation 



$(fc) 



3 

47r 



$(fc')dfc' 



(fc_/c')2+^2(fc2+j(,/2)2 



(46) 



In (02} we set a;"^ = a{T/uJo) with a w 3^/2^ to match Eq. 
H44II at vanishing r. To shorten notations, we dropped the 
subscript y from ky. 

Equation lli6|l can be analyzed analytically. Our goal is 
to understand why Tc drops by two orders of magnitude 
compared to Gwq already at relatively small r ~ 0.1. 

We search for the solution of II46|I in the form 
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""^-i^--^. 


: 


r 


3/4x - 1 

0.5 r-=^'= 


"5^ 


^ 


1 


" 



$(fc) = $o(fc)+'5$(fc), 



(47) 



where $o(fc) = fc is the unphysical solution at r ^ 
found in Subsection IIVAL and consider (5$(fc) as pertur- 
bation. Substituting into ijifiji and expanding in (5$, we 
obtain 



sm - II 



5<^{k')dk' 
(k - k'f + x^ 



3fc 



— -life 
4x 





3r2fc5 
2xq 



4a/ a;^ +47-2^4 



(48) 



where xq = 3/4, and we dropped the regular 0{r) terms 
from the integral 



<!>o{k')dk' 



(fc 



UA2 



^(F 



^2^1 2 



(49) 



and the higher-order terms in the expansion of the 
square- root in II48|I . We will see below that relevant k 
are of order r"^/^ and 3/4a; — 1 is of order r^/^, so that 



r2fc4 



^2/5 



which justifies both the expansion of the 
square- root in II48|I and dropping the regular 0{r) terms. 
We search for the solution of II48|I in the form 



J$(fc) = r^ {Ak^ + Bk^ 



\k\ < k„ 



(50) 



Figure 6: Scaling of the "critical temperature" x{r) (solid 
lines) and the position of the maximum fcm(f) (dashed lines) 
with the parameter r in the kernel 13611 of the one-dimensional 
"toy model". Numerical solution (red curves) closely follow 
the scaling laws obtained analytically (black straight lines): 
3/4a; -1^ 2r^/^ and fcm ~ r-^/^ 



Substituting this into lUsll and equating the prefactors for 
k, k^ and k^ terms in lUSJI . we obtain the set of equations 



Ax 



I A 



27r 



3 

ix 

(9fc„,A + 5klB) 



2xq 
15 



fcm-B 



TT 

3 

4a; 



Solving the set, we obtain 

(3/4a;-l)^ _ SOu^ ,^27 



(54) 
(55) 
(56) 

(57) 



where u — /cm(3/4a; — 1). A unique solution of this equa- 
tion exists when u = IB/tt and 



where fcm is a cutoff that has to be found from an aux- 
ihary condition that the solution (|17jl with 6^{k) in the 
form of l|5?11l is unique. Substituting lf5(Hl into lUsll and 
using 



4a; 



^2/5 



320*81 



1/5 



3^2/5 



(58) 



Using the definition of x we then obtain 



k'^dk' 



-K 



[k-k'Y + xl 

k'^dk' 
{k~k')'^ + xl 



we obtain 

3 

47r 
3r2 

47r 



S<P{k')dk' 
[k-k'Y+xl 

k{6k 



10 3 
A + -;rk^ 



Xo 



6kkn 



F— -|-20fc^fcn 
Xq 



B) 



k\-A 

Xq 



(51) 



yA;A;^(52) 



Tr-Q 



47r (1 -f- 3r2/5)4 



(59) 



20A;mB) 



k^^B 

Xq 



This formula is formally valid at small r, but it works 
surprisingly well for all r < 1 as evidenced by the com- 
parison of H59|l with the numerical solution of Ij48|l . We 
compare analytic and numerical results in Fig. El III 
Eq. IJ59II shows that Tc does indeed drop quite sub- 
(53)stantially, compared to the asymptotic linear behavior 
gH: at r ~ 0.1, (1 + ir'^/'^Y « 23. Furthermore, T^ 
given by ll59ll is rather flat at intermediate r ~ 0.04 — 0.5 
where Tc ^ 0.02g. This is consistent with the numerical 
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Figure 7: Solutions f{k) of the one-dimensional "toy model" 
for several values of the parameter r (r decreases for the curves 
from left to right). Dependence of the position of the maxi- 
mum kyn{r) is plotted in Fig. |^ 



solution of lUsll . The fiat beliavior in a wide range of r 
reproduces wiiat we iiave found numerically for the full 
gap equation (with the full frequency summation). 

We also see from ^^ that A ^ r~^/^, B -- r"^/^, and 
the cutoff fcm diverges as 



= ull 



\Ax 



1 



-2/5 



(60) 



term as the adjusting parameter. We found that the nu- 
merical data are reproduced reasonably well if we set this 
prefactor to 5 (instead of 3 in Eq. ll59ll ). We show the fit 
in Fig. 21 



C. Away from the QCP 



At deviations from the QCP towards larger dopings, 
i.e., into paramagnetic phase, Tc decreases and eventually 
disappears. The value of /cm, however, does not track the 
decrease of Tc, i.e., the d— wave gap extends over a finite 
momentum range along the FS even in the overdoped 
materials. At deviations into antiferromagnetic phase, 
the FS evolves into hole and electron pockets, and the 
locations of k^ gradually approach the locations of the 
hot spots. 

Our results for Tc and the gap survive even when the 
magnetic correlation length ^ remains finite at the doping 
where 2kp = Q and the Fermi surface has the form shown 
in Fig. 1. In this situation, the antiferromagnetic QCP 
shifts to lower dopings, when the hot spots are already 
away from the Brillouin zone diagonals. We found that 
the modifications of our results are small as long as Tc > 
J£,~^ where J ^ Vp/g is of the order of the the exchange 
interaction. 



Substituting fcm and A and B into Ij5n|l . we obtain 



S<^{k) ^ $o(fc) [{k/kn,f + a (k/k„,f) , a = Oil) 

(61) 
We see that at A; ~ /cm, the corrections to zero-order solu- 
tion become of order 1. A perturbation theory does not 
allow one to go beyond this scale (i.e., to analyze 6^{k) 
for k > fcm), but it is reasonable to assume that above fcm, 
<i)(fc) = $o(fc) + (5$(fc) begins decreasing. The numerical 
solution of lUsll confirms this (see Fig. CJ. We see from 
Fig. dthat the gap in the toy model is a non-monotonic 
function of kyi it is linear at small ky, passes through a 
maximum at ky ~ fcm, and then falls off. This behav- 
ior is again fully consistent with what we have obtained 
numerically for the full gap equation. 

Returning to the actual, unrescaled momenta ky, we 
obtain, using the definition of r, that the maximum of 
the gap is located at 



^1 

TTVp 



Aj.~2/5 
VF 



kpr 



3/5 



(62) 



where /3q — vp/kp. For /? ~ /3o, and r <ti I, k^ is much 
less than kp, i.e., the d— wave gap is confined to a near 
vicinity of the zone diagonal. 

We fitted the actual Tc for the full model by the func- 
tional form of Eq. Ij59|l , using the prefactor for the r^' ^ 



SUMMARY 



We considered the normal state properties and pairing 
near a 2kp antiferromagnetic QCP and applied the re- 
sults to electron-doped cuprates at optimal doping. We 
found that the breakdown of the Fermi-liquid description 
at QCP leads to peculiar frequency dependences of the 
conductivity and the i?2g Raman response. We found 
that Tc remains finite at the QCP. The pairing gap at 
QCP has dx2-y2 symmetry, but is highly anisotropic and 
confined to momenta near the zone diagonals. The value 
of Tc is about lOK for the same spin-fermion coupling 
as for hole-doped cuprates, where the pairing predomi- 
nantly involves antinodal fermions. The strong reduction 
of Tc compared to hole-doped case is due to the fact that 
for the d— wave pairing of near-nodal fermions, the cur- 
vature of the Fermi surface plays a major role. 

We acknowledge useful discussions with G. Blumberg, 
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Qazilbash, and V. Yakovenko. The research is supported 
by Condensed Matter Theory Center at UMD (P.K, A.C) 
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Appendix A: TRANSFORMATION OF THE 

POLARIZATION OPERATOR WHEN HOT 

SPOTS MERGE 

In this Appendix we detail the transformation of the 
polarization operator with the shift in the position of the 
hot spots along the Fermi surface. In hole-doped mate- 
rials the spots with the strong electron interaction are 
close to the (0, it) and symmetry- related points on the 
Fermi surface in the notations of Fig. Q] With doping 
and expansion of the Fermi surface its boundary crosses 
the antiferromagnetic Brillouin zone at points that be- 
come closer to each other and to the zone diagonal un- 
til in electron-doped cuprates they merge pairwise at the 
doping level when the Fermi surface just touches the anti- 
ferromagnetic Brillouin zone at {tt/2, it/2) (this situation 
is shown in Fig. QJ. In hole-doped materials it sufficed 
to take linear approximation of electron dispersion near 
the hot spots: ek — vpk. The polarization operator of 
non-interacting fermions then had the Landau damping 
form 



n(q,rj) = g'^\n\/'iTrv^u 



y> 



(Al) 



where g is the coupling, and vp — {vx,Vy)i^. In the sit- 
uation depicted in Fig. ^one has to keep second order 
terms in ek l|HJ, and the bare polarization operator has 
the form lllflll . Here we show how to go from IjAlll to 
ifTTHl with continuous change in the electron spectrum. 
For later use we will carry out slightly more general cal- 
culation keeping fermionic self-energy that depends only 
on frequency. Consider a spectrum that has both Hnear 
term in transverse momentum and curvature: 



Ek = v^kx + Vyky + 0^ky 
ek+Q = -Vxkx 



Vyky ■ 



P^kl. 



(A2) 
(A3) 



The case (3 — gives dispersion around well-separated 
hot-spots of hole-doped materials, and in the case Vy = Q 
we revert to IJHJ. To alleviate integration in k we note 
that the following equality holds 



Sk + Ck+q+Q 



where 



^q = 



-Vxqx 



P^q. 



2„2 

y 



2 2/32 

differs from Ec^ dll by -Vy/2(]'^, and 

^V2,(...|.^) 



(A4) 



- Td^ (A5) 



y 



(A6) 



is a convenient dummy variable in the integration instead 
of ky. The Jacobian of the substitution k ^ (ek, y) is 



9(ek,2/) 



c)[kx, ky) 



V2vxP. 



(A7) 



Momentum integration in (EJ then gives 



n(q,") = 257|^G(k, 



uj^'4)G(k + q,uj+'4) 



W 



47rV2wx/3 



sgn \uj -I- Y j + sgn \uj ~ ^ j 



y^q " * ( 



S , |fi| + S 



duj, (A8) 



where J^^i — ^ + S(w)- Simplifying frequency integration 
from the antisymmetric integrand 



ii{q,n) 



V2nvx(3 



du) 



^ 



n|/2 A/£^q +M5^,.,^m +^,.,_M 



UJ+^ ' w-J 



or, with the substitute u> ^ uj + \i^\/2, 



(A9) 



n(q,fi) 



9 



2'n:vxl3 Jo 



2 .00 WA/£;2 + (I]^+S^+|o|)2-^q 



i?2 + (I]^ + I]^+|^|)2 



-duj. 



(AlO) 

When T,{uj) = the integral above can be taken exactly 
and gives 



n(q,fi) 



n^+El+E^, (All) 



2'n:vx(i 

which coincides with IjlOII when Vy — and so i?q = 
Ec^. In the opposite Hmit /3 ^ 0, -Eq ^ -Vy/2f3'^ and, 
expanding ljAll|l up to first order in If^lZ-Bq gives 

I{{ci,n) ^ g^\Vt\/2TTVxVy, (A12) 

i.e. exactly twice the expression |IA1|I . accounting for the 
pair of hot spots that emerge. 

Appendix B: SELF-ENERGY AT T = 

In this Appendix we find self-consistent fermionic self- 
energy S(w) at low frequencies as explained in Section 
nil Bl Following the derivation that leads to Eq. Ijl0|l it 
can be shown that the polarization operator of interact- 
ing fermions that have self-energy S(ijj) is 



n(q,fi) 



9 



2ttvfP 



El + {^^ + ^^+\n\)''~E^ 



£2 + (I]^ + 2^+1^1 )2 



-duj. 



(Bl) 

where E^^ = w -f S(u;). It is easy to see that in the case 
E(a;) — the integral above can be taken exactly and 
gives ljl()|l . Note that we subtract divergent high-energy 
contribution 



Ho (0,0) 



9 



2ttvfP Jo 



duj 



2S, 



(B2) 
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because it is already included in the theory as the ^~^ 
"mass" term in the bare susceptibility IpJ- 

The above expression for 11 when substituted as the 
bosonic self-energy in the bosonic propagator ((Jjl in Eq. 
IjKill yields a self-consistent equation on E. We will seek 
a solution to this equation in the form of a power law: 
E^j ~ acu"', where cr is a constant. Using the substitutions 
il — \uj\t so that Ef2 = ^ojf^ , and uj' = \uj\ts, where 
uj' is the dummy variable in 11, so that Sj^' + Ecj'+n = 
So(s'^ + (s + 1)'^), and q = {^T^nf/'^v/ (5, we arrive at a 
condition on the exponent 7: 



TTiV 



1 



7 Jo 



f{v;-f)dv, 



(B3) 



where 



(B4) 



/(«;7) ~ Jo V V^4 + (57 + (s+ 1)7)2 

\/v/9w4 _^ (57 + (s + 1)7)2 _ 3^2 
V9w4 + (57 + (s + 1)7)2 

The right-hand side of Eq. IIB3|l is a continuous function 
of 7 in I < 7 < 1 . It monotonically grows from zero at 
7 = l^and crosses the constant nN/3^/2 at 7 ~ 0.85 (see 
Fig. |S||. Recalling that a used in the main body of the 
paper is a = 1 — 7, we get a « 0.15 as was ascertained. 
The calculation of the right-hand side of Eq. (IB3|I takes 
some care. Integral ljB4|l diverges at the upper limit when 
7- 



Expanding in large s gives for the residual of IJB4I 



.)ds 



v'sl-'-'/' 



-7/2 



^/2 \ 2(1 - 37/2) 



(B5) 



In 



whence J f{v;j)dv ^ y^T~~3jj2 for 7 close to 
the opposite Hmit 7 ^ 1 it is the integral over v in JB^ 
that diverges at the upper limit. At large v the scaling 
substitution s = v'^''^ z into IIB4ll separates an integral 
over z: 



1 



c'hy- 



so the residual of the integral over v m ^ 



fiv]l)dv: 

diverges as ~ (1 — 7)^^. The auxiliary integral 



(f-2)C(7) 



(B6) 



(B7) 



Cil) 



y^v/9 + (2^^)2 -'3 2 



x/9 + {2z-fy 



\dz. (B8) 




Figure 8: Graphical solution to Eq. llB3ll . Straight line and 
the curve represent the left and right-hand side of Eq. l|B3|l . 
Crossing is at about 7 ~ 0.85, or a = 1 — 7 ~ 0.15. 



It is negative in the interval 3 < 7 < 1 and diverges 
when 7 approaches |. When 7 = 1 the integral can be 
taken analytically: C(l) = — -s/G- To produce Fig. |Hlwe 
combined the numerical integration upto the finite cutoff 
with the analytical residuals found above. 



Appendix C: CHOICE OF THE CUTOFF FOR 
THE CONDUCTIVITY 



In this Appendix we substantiate our claim that the 
frequency dependence of the conductivity does not rely 
on how the integration over transverse momenta in Ij28|l 
is estimated. To evaluate the integral over ky we as- 
sumed that T,{e,ky) stays almost constant and equal to 



S(e) for |n,j/ 



K11 \ K' 



y,max 



and then falls off rapidly. Crudely 



we may just replace / dky with fcj,,max- However, since 
^y.max oc uj"' dcpcnds On frequency, it is not clear whether 
we should choose the external uj or the internal frequency 
e for the cutoff. Choosing the first one yields 



criiuj) ex Lo' 



7-1 



de 



Lu + S(w — e) -I- S(e) 



(CI) 



On the other hand, if we simply integrate over ky in Ij28|l 
with S(e, /cy) = S(e) for \ky\ < fcy,max, the result is 



anuj} 



ex 2 



V2 





{CO 



-T.{lu 



E(e) 



w + S(a; — e) 



de. 



(C2) 



Intuitively, the choice of the cutoff of the momentum in- 
tegration should not alter the frequency dependence of 
the conductivity since the internal e stays of order of the 
external lu anyway. And indeed we found that all three 
estimates lead to the same pseudo-scaling of cr(w) in the 
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range ljq < to < 40wo- The choice of the cutoff alters 
only the irrelevant preexponent shared by both "^a and 

5(7. 

The fact that ^a and '^a scale together with |cr((jj)| 
can be seen from arctan(3cr/5Rcr) staying almost con- 
stant over the range ujq < oj < 40wo (see Fig. [2l right 
panel). The value of the constant f = arctan(3(7/3ficr) 
itself changes with the Wick transition from Matsubara 
to real frequencies. In particular, ip depends on whether 
the Matsubara or real frequency enters the momentum 
cutoff A:y^max(w). We argue that it is logical to cut off the 
momentum integration on Matsubara frequencies since 
the subsequent frequency integration is done over Mat- 
subara frequencies as well. With this choice the constant 
(/? « 1 (see Fig. [21 right panel). 



Introducing S„ 
tion on A„: 



: OJnZn, ^r 



A„Z„ we get one equa- 



A„ = T V" I — ^ I SgTLUJmXi^n - Wm) (D7) 

^ V ^m W„ ) 

Here the summation clearly can be restricted to m ^ n. 
Returning back to the original formulas we see that we 
can redefine S„, $„ as sums with m ^ n. (For the self- 
energy in the form ljD8|l this means ri„ ^ 0). 

Expressing Matsubara frequencies in overt form, 
bosonic as Vtn = 2TTTn, fermionic as a;„ — 2nT{n+ 1/2), 
we may rewrite the finite sum in the self-energy IJD3II via 
the harmonic numbers 



Appendix D: ELIASHBERG LOCAL PAIRING 

In this Appendix we describe the numerical solution 
of the pairing problem II42II for a general power-law local 
interaction. The results of this Section were used for Tc 
of (EHJ. 

The general interaction is (the frequencies and the tem- 
perature are measured in suitable units which we put 
unity) 



x(f^„) = ^(i-7)ia.r 



(Dl) 



Fermionic self-energy S(w„) corresponding to ljDl|l is 
given by 



S(tJ„) = r^SgnLJ„x(^n - ^m). 



(D2) 



On substituting fin = ujn — uJm the sum in S„ becomes 
finite: 

Y.{uJn) = sgnuJnT J2 Xi^m)- (D3) 

|0„|<|cj„| 

Approximating the sum by an integral we get the limiting 
scaling of the self-energy: 



1 r 
E(w) « sgnw— / x{^)d^ 
71" Jo 



sgnw|a; 



1-7 



(D4) 



Reverting to the pairing problem and denoting S]„ = 
u!n + S(w„) and $„ = $(w„), we get for the linearized 
system on the self-energy and anomalous vertex at finite 
temperatures by analogy to Eqs. II39|I . II43|I : 



-T^SgnUJmXi^n-l^Tn) (D5) 



$n 



^E 



$. 



V 



■X(Wn - UJm) 



(D6) 



1 

i/(7) = yJ- 

m— 1 



(D8) 



as 



E(c^„>0) = (l-7)(2^r)i-^if(^). 



(D9) 



Eq. on the anomalous vertex $ should be under- 
stood in the sense that at some T^ which would be the 
critical temperature for the onset of pairing, a non-trivial 
'^{iOn) would become a solution, i.e. at T^ the largest 
eigenvalue of () crosses unity. The kernel of Eq. () 
symmetric with respect to simultaneous change of sign 
before both w„ and w^, so its solutions must be ei- 
ther even or odd in frequency. From a known theorem, 
the largest eigenvalue corresponds to an even eigenvec- 
tor <J>(w„) = $(— cj„). For the even solution we get the 
equation 



$n = 



E' 



TO 



■ TO 



(2^T) 



-(TO + i) + ir, 



(7) 



■$„ 



(DIO) 



where in the term m ^ n expression \n- 
substituted with zero. 



-toI '>' should be 



The numerical solution of the eigenvalue problem was 
done with LAPACK for a finite matrix of order TV with 
several increasing values of N. The sequence of the criti- 
cal temperatures Tc{N) found for each N is then extrap- 
olated to 1/A^ ^ 0. A sequence for 7 = 1/4 plotted in 
Fig- ©) right panel, represents a typical picture. The 
critical temperature as function of 7 are shown in left 
panel of Fig. EI}. 




power law exponent, 7 



0.02 0.04 0.06 0.08 
inverse matrix order, 1/N 



Figure 9: Critical temperatures Tc found numerically vs. the 
power law exponent 7 (left panel). Right panel represents 
the dependence of Ta on the matrix size A'^ extrapolated to 
1/A^ — > for 7 = 1/4; behavior for other 7 is similar; left 
panel shows Tc already extrapolated to 1/N -^ 0.} 
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